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ABSTRACT 

We consider we problem of determining of D-optimal experimental designs. During many possible experimental 
plans we consider designs called chemical balance weighing designs that are described by the linear model. We study the 
problem of determining optimal design under assumption that the measurement errors are equal correlated. 
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1. INTRODUCTION 

The planning of experiments is very important part of mathematical statistics. Among many conditions 
determining best experimental plan, very significant is to fix the optimality criterion. The choice of optimality criterion 
influences on the properties of determined estimators. In many papers, some problems related to the determining of the 
optimality criterion and the consequences of this choice are considered (See Pukelsheim, 1993). For the most part, the 
optimality criteria there are functions of the design matrix. In determining the best experimental plan, we have to 
remember about the experimental model. In present paper we consider chemical balance weighing design. Such designs 
can be applied in economy, agriculture, optics. Chemical balance weighing design there is the model of experimental 
design that permits to determining unknown measurements of p objects if we will take them n times in different 

combinations to measure. The results of n measurement operations of p objects will fit into the linear model 
y = Xw + e, where y is the nxl random vector of observed weights, Xe$ ropm (-l,0,l), O nx m (— 1, 0, l) 
denotes the class of nx p matrices X = [XyJ, i = 1,2,.../! , j = 1,2,..., p , of known elements equal to —1, 1 or 
according as in the i th weighing operation the j th object is placed on the right pan, left pan or not and 
m = maxjmj, m 2 , ...,m p J, where trij, j ; = 1,2, ...,p, represents the number of elements equal to —1 and 1 in j th 
column of X , W is a p x 1 vector representing unknown weights of objects and e is an nxl random vector of errors. 
Additionally, we assume that there are not systematic errors, i.e. E(e) = n and Var(e)=<T 2 G, where (J 1 is the 
constant variance errors, G is nxn known positive definite matrix of the form 

G = g[(l-p)l„+pl,X] g>0, — <P<l (l) 

n — l 
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The form of the matrix G G means that the measurement errors are correlated. Let note, 



G = 



I L i i 

n l + p(n-l) n n 



For estimation W we use the normal equations XG 'Xw = XG l y , 



Any chemical balance weighing design is singular or nonsingular depending on whether XG 'X is singular or 
nonsingular, respectively. Since G is known positive definite matrix then X G X is nonsingular if and only if X is of 
full column rank. However, if XG *X is nonsingular then the generalized least squares estimator of W is given by 

W = (X'G _1 X) _1 X'G -I y and Var(w)=cr 2 (x'G _1 x) \ The matrix M = X G 'X is called the information 
matrix of the design X . 

In many problems concerning experimental designs the D-optimal designs are considered. In the case of chemical 
balance weighing designs, the design X rf is D-optimal in the class of the design matrices *F (Z O nxp m (— 1, 0, l) if 

det(X rf G _1 X rf )= max{det(M): XeT}. It is known that det(M) is maximal if and only if det(M _1 ) is 

minimal. The concept of D-optimality was considered in the books of Raghavarao (1971), Banerjee (1975), Shah and 
Sinha (1989). In the paper Jacroux et al. (1983) the problems related to the D-optimality were presented for the case 

G = \ n . In Masaro and Wong (2008), Katulska and Smaga (2013), D-optimal weighing designs X e A nxp (— 1, l), 
where A nxp (— 1, l) is the set of all nx p matrices with elements —1 or 1, only, are considered. 

In this paper, we present new results concerned D-optimal chemical balance weighing designs that are the 
generalization of given in Katulska and Smaga (2013) ones. We assume that the random errors are equally correlated and 
they have equal variances. We give lower bound for the determinant of the inverse of information matrix. Also we give the 
examples of the D-optimal chemical balance weighing design for which the determinant of the inverse of information 
matrix attains the lower bound. 

2. D-OPTIMAL DESIGN 

Let X = [x l5 X 2 , ...,X ] G$ (-1, 0, l) . From Section l.c.l(ii) (b), Rao (1973) we have 

Lemma 2.1: Let G be as in (1). For diagonal elements of the inverse of information matrix the following 
inequality holds 

M 'x/) ' > ^ ■ GO 

XjXj \ + p{n-\) 



Now, we prove the inequality which gives the lower bound for the determinant of the inverse of information 

matrix. 



Theorem 2.1: Let G be as in (1). If XeO m (-l,0,l) then 
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det(M _1 ) 



> 1 



V m y 
8(1-P) 



if 0<p<l 



m — 



p{m-2uf 



(3) 



if 



n-1 



<p<0 



\ + p(n-l) 

Where u =1X11111^^ U 2 , j, U j represents the number of elements equal to —1 in j th column of X, 

y' = l,2,...,/?. 

Proof: By the Hadamard's inequality it follows that 

det(M- 1 )>nM]] = n(x' 7 .G- 1 x,)r 1 . (4) 

7=1 7=1 

Hence by Lemma 2.1 we obtain \XjG~ l Xj] > g(l — p^XjXjj taking Xjl n = 0. The maximal value of 
X -X ■ is m assuming that X e <& n> , pm (—l, 0, l). For the case — (ft — l) 1 <p< , the maximal value of denominator 

X i^/i^ / P ' ' i \2 

is attained if -, r- is maximal, so if X -l n l n x • = [m — 2u). So, the proof is completed. 

l + p(n-l) 

Definition 2.1: Any chemical balance weighing design X G <& nxp m (— 1, 0, l) with the covariance matrix of 
errors <J 2 G , where G is given by (1), is said to be regular D-optimal if it satisfies the equality in (3), that is 



det(M -1 )=< 



g(i-pY 



y m 



if 0<p<l 



r(l-p) 



m- 



p(m-2uf 



if 



1 



(5) 



n-1 



<p<0 



l + p(n-l) 

Let note, the regular D-optimal design is D-optimal, whereas the inverse sentence may not be true. 
Theorem 2.2: Any chemical balance weighing X G <D nxp m (— 1, 0, l) with the covariance matrix of errors 
(7 G , where G is given by (1), is regular D-optimal if and only if 
XX = mI p if p = 0, 

X X = ml p and X'l„ = p if < p < 1, 



XX = mI p - ; V I p -l p l p andXl^z p if -<p<0, 

1 + p[n - 1) n-l 
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Where Z is pxl vector for which the j th element is equal to (m — 2u) or — (m — 2u) , j = 1,2, ...,/?. 
Proof: Let X e *& nxp m (— 1, 0, l) be the design with the covariance matrix cr 2 G , where G is given by (1). 

• If XX = mI p then X -X^ =0 for j ^ k , j, k =1,2,...,/? . Because of this and from Lemma 2.1 it follows 

that if p = then the inverse of the information matrix is equal to — I , so its determinant satisfies (5). 

m 

• If < p< 1 and x 7 -l n =0, 7=1,2,...,/?, then under Lemma 2.1 we deduce that the inverse of the 
information matrix of the design X is equal to g(l — p)m [ l p and therefore the equality (5) is also true. 

■ p{m-2xif ( < \ 

• If XX = mI -— — 7 M I -1„1„ and XI =Z„, then by Lemma 2.1 we deduce that the inverse 

p \ + p{n-\V p p p > n p 



v U i p(m-2u) 
of the information matrix of the design A is equal to g\Y — p\ m — -, 

\ + p\n-\) 



2 V 1 



\ p and therefore the 



equality (5) holds. 

• Thus the design X is regular D-optimal in many cases. Hence Theorem. 

Corollary 2.1: Let < p < 1. The necessary condition for existence regular D-optimal chemical balance weighing 
design X e *f> nxp m (— 1, 0, l) with the covariance matrix of errors <T 2 G , where G is given by (1), is m = 0mod(2) . 

3. EXAMPLES OF REGULAR D-OPTIMAL DESIGN 

In order to determine regular D-optimal chemical balance weighing design XeO nx? m (-l, 0, l) with the 

covariance matrix of errors (7 G , where G is given by (1), we have to fix regular D-optimal design in the class 

m (— 1, 0, l) for each p separately. It means for p = 0, < p < 1 and for < p < . The assignation of 

n-l 

the design matrix that satisfies one of the conditions 1-3 from Theorem 2.2 is connected with determining the matrix with 
special properties. Some methods of construction of such matrices given in literature is based on the incidence matrices of 
some known block designs, for example balanced incomplete block designs, ternary balanced block designs, balanced 
bipartite block designs. In different papers there are given the algorithms nominated the best design (See Angelis et al., 
2001). The aim of this paper is to give the lower bound of the determinant of the inverse of information matrix for the 
design and determining of the conditions under that the lower bound is satisfied. So we will give only some examples of 
regular D-optimal designs. 

Let p = 0. We consider X e €> 10x5 8 (-l, 0, l) . We have Var(w, )= — (J 2 , J=l,2,...£. Regular D- 

8 

optimal chemical balance weighing design in this class is given in the form 
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Now, let us 


consider 


XeO> 
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weighing design for any p , < p < 1 and Var (\Vj ) = 




- p K 
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< p < we consider G = g 
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. So, 



. So, we determine the regular D-optimal 



design in the class X e ^loxs gl - 1, 0, l) . We obtain 
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And Var(w ; .)=^j-a 2 , j = l2,...$. 
4. DISCUSSIONS 

We consider the regular D-optimal chemical balance weighing design in order to determine estimators of 
unknown measurements of objects with the smallest variance. So, let us consider X e *& nxp m {— 1, 0, l) for fixed number 



of n , p , m . We consider parameter p in intervals given in Theorem 2.2. If p = then Varfw ■ )= — CT 2 , if 

m 
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Q<p<\ then Var(w)=^^a 2 if ^-<p<0 then Var(w,)=g(l-p{m- ^ m 7 2M)2 J 

m n-l J/ I l + p(n-lj 



cr 2 



7 = 1,2,...,/?. 

In the planning of each experiment according to the model of the chemical balance weighing design, as the first 
step we have to choice the value of correlation of errors, i.e. p . Let us consider any p t , < p < 1 , t = 1,2 , p l ^ /? 2 • 

It is worth pointing out that the design X G *& nX p,m i~ 1> 0' l) satisfying Theorem 2.2 is regular D-optimal in the sense of 

attaining equality (5), for and p 2 . Thus the design Xe$ nx;) m (-1, 0, l) is regular D-optimal for any /? , 

< p < 1 . Simultaneously the lower bound of variance of estimators is not the same for different numbers of p . 
According to the paper Masaro and Wong (2008) the design is robust. 

The choice of the value of p is very important from the point of view of estimation of unknown measurements of objects. 
For < p < 1 , greater value of p implies minor value of Var (w • ) . For — (ft — l) 1 < p< , minor value of p 
implies major value of Var {\Vj ) . 

In the case when the regular D-optimal design doesn't exist we have to determine the design for which the 
determinant of the inverse of the information matrix is the closest to given in (5) one. The lowest bound given in (5) 
permits to check if given design is optimal or not. In such cases some problems of determining optimal designs are given in 
Graczyk (201 1,2012). 
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